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1  a  ABSTRACTS*"*^  W«®d< 

This  report  describes  a  track-co-track  fusion  alj^orithm  for  airborne  surveillance 
systems  employing  multiple  dissimilar  sensors  (radar,  IR,  and  laser  radar).  I'hese 
sensors  detect  targets  and  create  tracks  at  different  data  rates.  The  algorithm 
presented  performs  synchronization  by  predicting  the  slower  tracks  to  the  update  times 
of  the  faster  tracks.  The  synchronized  tracks  are  then  tested  for  association  to 
determine  whether  or  not  the  two  tracks  originated  from  the  same  target.  It  is  shown 
that  the  probability  distribution  of  correct  track  association  can  be  improved  if  the 
test  statistic  for  association  incorporates  cross-covariance  between  the  two  tracks. 

A  recursive  algorithm  for  computing  the  cross-covariance  is  obtained.  In  addition,  a 
trade-off  study  involving  the  probability  of  correct  association,  number  of  track 
matching  points,  size  of  association  gate,  and  probability  of  false  correlation  has 
been  prepared.  These  algorithms  are  coded  in  MATLAB  and  the  results  of  simulations 
confirming  the  proof  of  concept  are  also  presented. 
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SECTION  1 


INTRODUCTION 

Design  of  current  tactical  air  defense  systems  rely  on  inputs  from  a  network  of  distributed  sensors 
and/or  platforms  to  construct  an  air  picture  and  to  coordinate  vital  mission  activities  such  as 
surveillance,  early  warning,  and  weapons  control.  Integration  of  various  data  channels  enhance 
effectiveness  of  surveillance,  improve  immunity  to  interference  and  deception,  increase  information 
capacity,  and  reduce  eclipsed  (where  data  are  unavailable  because  of  eclipsing  by  the  transmission) 
and  blind  zones.  Current  system  architectures,  however,  consist  of  separate  track  files  being 
maintained  by  individual  sensor  platforms.  This  track  data  must  be  associated  and  fused  to  provide  a 
coordinated  assessment  of  the  tactical  environment.  The  de  facto  standard  for  the  exchange  of  air 
surveillance  data  is  through  the  use  of  track  files  over  non-dedicated  tactical  data  links.  The  contents 
of  these  messages  limit  the  data  available  for  track-to-track  processing,  and  incompatible  network 
protocols  constrain  the  dispersion  of  track  files  among  users.  As  a  result,  existing  tactical 
surveillance  systems  are  unable  to  provide  an  integrated  air  picture  obtainable  in  a  centralized  fusion 
system  architecture. 

Previous  work  on  track-to-track  fusion  was  based  on  the  assumption  that  the  estimation  errors  of 
tracks  for  the  same  target  obtained  from  different  sources  are  uncorrelated  [1].  This  assumption  is 
incorrect,  however,  because  the  process  noise  associated  with  target  maneuvers  is  common  to  the 
filter  dynamics  used  by  the  estimators  [2].  Hence,  the  target  tracks  are  correlated  even  though  the 
sensor  measurement  errors  are  uncorrelated.  The  effect  of  the  cross-covariance  on  the  performance  of 
track-to-track  association  and  kinematic  fusion  was  investigated  in  [3]  for  fusion  of  tracks  created  by 
two  dissimilar  sensors.  This  problem  was  also  investigated  in  [4]  for  fusion  of  synchronous  tracks 
from  two  identical  sensors.  However,  results  published  in  unclassified  literature  do  not  discuss  fusion 
of  asynchronous  tracks. 
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The  goal  in  this  contract  is  to  fuse  asynchronous  tracks  which  are  updated  at  different  time  update 
rates.  Specifically,  a  Kinematic  Track  Synchronization  (KTS)  algorithm  for  synchronizing  fast  and 
slow  tracks  was  developed.  Then  the  minimum  closeness  score  criteria  was  used  in  the  Kinematic 
Track  Association  (KTA)  algorithm  to  test  the  hypothesis  that  two  candidate  tracks  for  fusion 
originated  from  the  same  target.  The  associated  track  pair  is  then  combined  using  a  Kinematic  Track 
Fusion  (KTF)  algorithm 

l.l  BACKGROUND 

Most  currently  fielded  surveillance  systems  employ  ad  hoc  methods  for  correlating  tracks  from 
multiple  sources.  These  methods  are  dependent  upon  specific  system  configurations  and  often  lack 
rigorous  theoretical  foundations.  Typically,  little  is  done  to  fuse  tracks  which  originate  from  the  same 
target.  Rather,  a  preferred  source  of  track  data  is  usually  identified,  and  the  tracks  from  that  source 
are  used  to  represent  the  correlated  aggregates,  while  the  rest  of  the  track  file  is  ignored.  A  track 
quality  (TQ)  parameter  provided  by  each  track  file  is  often  the  discriminant  used  for  selecting  a 
preferred  track  source.  Unfortunately,  there  is  no  rigorous  criteria  for  determining  TQ,  instead,  all 
systems  are  free  to  compute  TQ  as  desired.  This  processing  approach  is  adopted  because  of 
constraints  imposed  by  the  requirement  to  adhere  to  existing  tactical  data  link  protocols 

The  Air  Force  Scientific  Advisory  Board  Summer  Panel  for  1991  recommended  a  theater  wide 
information  exchange  architecture  for  distributing  tactical  surveillance  data  in  track  form.  To 
accomplish  this  task,  a  new  breed  of  real  time  RF  links  is  being  planned  to  link  a  vast  array  of 
platforms,  e.g.,  E-3  AWACS,  E-2C  Hawkeye,  RC-135  Rivet  Joint,  Common  Sensor  Platforms  (RC- 
12  Guardrail,  EP-3  and  ES-3A),  E-8C  Joint  STARS,  FLIR/TV-equipped  UAVs,  etc.  This  new  breed 
ot  battlefield  network  links  sensors  and  weapons  to  make  it  possible  not  only  to  launch  missiles  or 
vector  aircraft  against  targets  far  beyond  the  range  of  each  weapon's  own  guidance  systems,  but  also 
enables  sensors/  platforms  to  share  information  and  automatically  update  the  information  in  a 
multitude  of  networked  fire  control  computers.  Put  simply,  this  kind  of  real-time  C^I  systems  not 
only  could  coordinate  sensors  and  shooters,  but  they  also  could  provide  a  rich  easily  accessible  body 
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of  tactical  information,  updated  in  real  time,  to  give  every  command  center  a  coherent  situational 
picture  of  a  conflict. 

1,2  APPROACH 

This  study  pursues  the  development  of  algorithms  which  can  integrate  all  available  track  information 
into  composite  tracks  of  greater  quality  than  the  constituent  parts.  The  principle  thrust  of  this  contract 
is  the  derivation  of  track  synchronization,  track-to-track  association,  and  track  fusion  algorithms. 

In  a  surveillance  system,  targets  are  detected  by  different  sensors  which  exploit  different 
characteristics  of  the  optical,  infrared,  and  electromagnetic  spectrum.  After  detection,  report  to  track 
association  is  performed  and  a  tracker  is  employed  to  create  tracks.  These  tracks  are  created  at 
different  rates  which  depend  upon  the  sensor  characteristics.  When  these  asynchronous  tracks  arrive 
at  the  fusion  center,  there  is  a  need  to  synchronize  them  prior  to  fusion.  In  this  report.  Kinematic 
Track  Synchronization  (KTS)  is  performed  by  predicting  the  slower  tracks  to  the  time  of  arrival  of 
faster  tracks.  In  addition,  covariance,  of  the  predicted  track  is  also  computed  which  is  used  in  the 
track  association  and  fusion  algorithm. 

A  difficult  aspect  of  track  level  fusion  is  the  decision  process  whereby  tracks  from  different  sources 
are  determined  to  represent  the  same  target.  In  this  report.  Kinematic  Track  Association  (KTA)  is 
performed  by  employing  a  track  matching  algorithm.  This  algorithm  employs  an  association  gate 
around  faster  tracks  and  performs  a  chi-square  test.  Slower  tracks  passing  this  test  are  then  used  to 
compute  closeness  score.  The  track  pair  which  has  the  minimum  closeness  score  is  then  selected  for 
association. 

In  this  report,  it  is  assumed  that  all  remote  tracks  are  transmitted  to  a  centralized  fusion  center.  At  the 
fusion  center.  Kinematic  Track  Fusion  (KTF)  is  performed  and  the  central  track  library  is  updated  by 
the  newly  created  fused  track.  In  this  report,  track  fusion  is  performed  by  employing  a  weighted 
covariance  algorithm  [3].  For  two-sensor  track  fusion,  this  algorithm  combines  associated  tracks  by 
weighting  each  of  these  tracks  by  certain  covariances  which  are  fimctions  of  their  respective  filtering 
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variances  as  well  as  their  cross-covariance.  An  algorithm  for  computing  the  cross-covariance 
between  asynchronous  tracks  is  also  derived. 

The  KTS,  KTA,  and  KTF  algorithms  developed  in  this  contract  are  implemented  in  MATLAB. 

These  algorithms  are  tested  by  creating  an  Airborne  Surveillance  Scenario  which  consists  of  four 
targets  which  are  flying  at  a  speed  of  0.3 1  km/sec.  maintaining  a  constant  altitude  of  9. 1  km.  These 
targets  are  sensed  by  three  sensors  (Radar,  IR,  and  Laser  Radar)  which  are  colocated  and  stationary  at 
9. 1  km.  Sensors  and  the  scenario  are  simulated  using  MATLAB.  A  Graphical  Lfser  Interface  (GUI) 
is  also  designed  in  Visual  Basic. 

1.3  OVERVIEW 

This  report  is  organized  as  follows.  In  section  2  mathematical  model  for  track  fusion  is  discussed. 

The  KTS  algorithm  is  described  in  section  3.  In  section  4  details  of  the  KTA  algorithm  is  given. 
Closed-form  mathematical  expressions  for  the  probabilities  of  correct  association  and  false 
correlation  are  derived.  In  addition,  trade-off  studies  involving  several  system  parameters  are  also 
presented.  Description  of  the  KTF  algorithm  is  given  in  section  5.  This  algorithm  incorporates  cross¬ 
covariance  between  asynchronous  tracks  which  originated  from  the  same  target.  In  section  6  details 
of  numerical  simulation  of  the  fusion  algorithm  for  an  airborne  surveillance  scenario  are  presented. 
Finally,  summary  of  the  project  and  some  concluding  remarks  are  contained  in  section  7. 


4 


SECTION  2 


MATHEMATICAL  MODEL 

For  the  sake  of  simplicity,  it  is  assumed  that  two  sensors  of  dissimilar  quality  (characterized  by 
differences  in  their  measurement  variances)  are  tracking  the  same  target.  After  acquisition,  each 
target  is  tracked  with  a  Kalman  filter  associated  with  the  sensor.  The  mathematical  model  describing 
the  target  dynamics  is  assumed  to  be  linear  and  of  the  form 
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(2-1) 


=<t)X(tk)  +  GW(tk) 


where,  the  target  state  vector  can  be  expressed  by  a  six  dimensional  vector 
X  =  (u^,  v^]^  with  the  three  dimensional  position  vector  u  and  three  dimensional  velocity  vector  v. 
The  input  noise  W(t|{)  is  assumed  to  be  a  three  dimensional  zero  mean  Gaussian  noise  vector  whose 
variance  is  Q  =  ql  where,  I  is  a  3x3  identity  matrix  and 


GE[W(tk)W(tO^]G^=C}QG^  = 


AT^  ,  AT^ 


3 

AT^ 


ql  ^qi 

ql  ATql 


(2-2) 


where,  E[.]  denotes  statistical  expectation.  This  target  is  tracked  by  three  sensors  whose  measurement 
model  is  described  by 


Z‘(tk)  =  h‘(X(tk))  +  v‘(tk) 


(2-3) 


5 


where,  the  superscript  i  denotes  the  type  of  sensor.  In  this  report,  three  different  sensors  (Radar,  IR, 
and  Laser  Radar)  are  used  to  track  the  target.  The  observation  vector  Z  (tiJ  for  Radar  is  [R,  A,  E],  for 
IR  is  [A,  E]  and  for  Laser  Radar  is  [R,  R,  A,  E]  where,  A  denotes  Azimuth,  E  denotes  Elevation,  R 
denotes  Range,  and  R  denotes  Range  Rate.  First  partial  derivatives  of  the  observation  matrix 
h(X(tk)),  are  given  below: 

For  radar  sensor. 


cosAcosE  sinAcosE  sinE  0  0  0 
-sinA  /  RcosE  cosA  /  RcosE  0  0  0  0 
-cosAsinE  /  R  -sinAsinE  /  R  cosE  /  R  0  0  0 


where,  (R,  A,  E)  is  the  predicted  measurement  vector. 
For  IR  sensor, 


-sinA  /  RcosE  cosA  /  RcosE  0  0  0  0 
-cosAsinE  /  R  -sinAsinE  /  R  cosE  /  R  0  0  0_ 


where,  (A,  E)  is  the  predicted  measurement  vector. 
For  Laser  Radar, 


cosAcosE 
-RcosAsinE  /  R 
-sinA  /  RcosE 
-cosAsinE  /  R 


sinAcosE 
-RsinAsinE  /  R 

cosA  /  RcosE 
-sinAsinE  /  R 


sinE  0  0  0 

Rsinl/R  cosAcosE  sinAcosE  sinE 

0  0  0  0 

cosE  /  R  0  0  0 


where,  ( R,  R,  A,  E)  is  the  predicted  measurement  vector.  The  measurement  noise  vector  for  the  three 
sensors  are  assumed  to  be  zero  mean  white  Gaussian  noise,  whose  variances  are 
described  by 
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R‘^  =  diag((T|,a^^,ff|)  ,  r' =  diag(c^,a|)  ,  R‘'  =  diag((a|,a|,aA,<y|) 


The  state  estimates,  filter  gains  and  covariances  of  the  tracks  are  obtained  from  the  following 
Extended  Kalman  Filter  (EKF)  equations: 

X'"(tk.,/tk.,)  =  <i»X'"(tk/tO  +  K'"(tk,,)[Z'"(tk,,)-H'"(t,^,)<|)X'"(tk^  (2-4) 

K"  (tk+i)  =  P”  (tk+i/t|c)  H”  (I»H)P"'(tk+l/tk)  (OxVl-'  (2-S) 

P'"('k+l/'k)*  +  GQGT  (2-6) 

P"  Ok+l/tk+1)  -  ( I  -  K”  (tk+1)  H"  (v,)l  P”  (<k+l/«k)  P-7) 

where,  the  superscript  m  denotes  the  type  of  sensor  which  created  the  track  and  it  is  assumed  that  the 
filter  dynamics  are  the  same  as  that  of  the  target.  It  is  also  assumed  that  these  tracks  are  created  at  the 
remote  stations  every  AT  sec.  and  transmitted  to  the  fusion  center  over  a  data  link  at  the  same  rate. 
Time  lags  and  delays  are  neglected.  The  input  noise  W(t,J  introduces  cross-correlation  between  the 
trackers  which  are  tracking  the  same  target.  Effect  of  the  cross-correlation  between  synchronous 
tracks  is  reflected  in  the  cross-covariance  matrix,  denoted  by  P' ,  which  is  given  [2]  by 

P<:(tk/tk)-E[X(tO-  X'(lk))[X(y-  Xi(tk)f 

=  [  I  -  K'  (tk)  H'  H  P'Ctk-l/'k-l)  [  I  -  K'(«k)  H'lT 

+  [  I  -  K'  dk)  H'  1  GQGT  1 1  •  K*  (tk)  tflT  (2-8) 

with  P<^(0/0)  =  0.  An  estimate'  of  the  fused  track  and  its  covariance,  respectively,  is  given  by  [3] 

X‘’(ti/tk)  =  X'(tk/tk)+[P'(tk/tk)-P'=(tk/tk)l(P')-'[X'(tk/tk)-X'(tk/tk)l{2-9) 

'  Recently,  it  has  been  shown  [6]  that  this  estimate  is  not  optimal 
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where 

P^(tk  /  tk )  =  P‘(tk  /  tk)  +  P^(tk  /  tk)  -  P'(tk  /  tk )  -  P'Ok  !  tk  f  (2-1 0 

At  steady  state,  equation  (2-8)  can  be  written  as 

pc  =  Fi  pc  F2^  +  Q* 

where, 

Fi  =  [  I  -  K‘h‘  ]  <j) 

F2*[I*KjH'](|>  (2-13) 

and  the  Kalman  filter  gains  for  the  sensors  are  described  by 

(2-14) 

In  addition,  the  input  matrix  Q*,  shown  in  (2-12)  is  given  by 

Q*  =  ( I  -  K'  H'  ]  GQGT  [  I  -  K  HMI"  (2-15) 

Equation  (2-12)  is  an  asymmetric  discrete  Lyapunov  equation  whose  steady  state  solution  was 
obtained  by  means  of  matrix  inversion  in  [5]. 
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SECTION  3 


KINEMATIC  TRACK  SYNCHRONIZATION 

In  today's  Command,  Control,  Communication,  and  Intelligence  (C^I)  environment,  remote  stations 
are  required  to  transmit  their  local  tracks  instead  of  raw  measurements,  to  a  fusion  center  over  a 
communication  link.  The  communication  links  have  different  chaimel  capacity,  bandwidth,  and  data 
rates.  Hence  the  remote  tracks  arrive  at  the  fusion  center  at  different  rates.  The  track  fusion 
algorithm  implemented  in  this  contract,  takes  into  account  the  rates  of  transmission  for  each  link. 
Specifically,  slower  tracks  are  predicted  to  the  arrival  time  of  the  faster  tracks.  In  addition, 
covariance  of  the  slower  track  is  updated  by  propagating  it  as  follows; 

X‘(tk.m/tk)='t>'"X‘(tk/tk) 

and  (3-1) 

P‘ (t|„„  /  )=♦  "P* (tk  / 1|, )»“')%!;♦ '■'OQG^tt  j-'r 


where,  the  superscript  i  denotes  the  track  number,  t|c+,„ *  (k+m)AT,  and  (jj™  =  It  is  assumed  that 

the  track  i  was  last  updated  at  time  .  When  this  track  is  associated  with  another  track)  at  time  t^+m , 
the  association  algorithm  takes  into  account  the  updated  covariance  of  the  propagated  tracks  as  shown 
in  (3-1).  Detail  of  the  association  algorithm  is  described  in  the  next  section.  If  the  association 
algorithm  accepts  the  hypothesis  that  the  two  tracks  originated  from  the  same  target,  then  the  track 
fusion  algorithm  combines  the  track)  and  the  propagated  track  i.  Detail  of  the  track  fusion  algorithm 
are  given  in  section  5. 
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SECTION  4 


KINEMATIC  TRACK  ASSOCIATION 

The  Kinematic  Track  Association  (KTA)  algorithm  implemented  in  this  contract  involves  testing  the 
hypothesis  (at  every  time  update  t)()  that  the  two  tracks  X'  and  originated  from  the  same  target  [7], 
Hence,  the  null  hypothesis  can  be  stated  as:  Hq:  X'  —  X^  =0  vs.  Hj:  X'  -X^  ;*0.  If  the  tracks  are 

independent  [I  ],  then  under  hypothesis  Ho,  cov(X‘  -  X^ )  =  PU  Assuming  Gaussian  probability 
distribution  of  the  track  estimates,  the  test  of  hypothesis  can  be  restated  as: 

Ho:(X‘  -Xj)'^(P‘  +pj)"‘(X‘  -Xj)<Xn  (4-1) 

where,  denotes  the  gate  width  of  the  association  region.  Since  the  test  statistic 

(X'  -  X  j  )^(P‘  +  P^  )■*  (X'  -  X()  is  chi-square  with  n  degrees  of  freedom,  the  test  threshold  X„  can 

be  chosen  for  a  level  of  significance  a  such  that 

Pr{(X‘  -X0^(P‘  +P^)“'(X‘  -XJ)>A.„  /Ho}=a  (4-2) 

4.1  PROBABILITY  OF  CORRECT  ASSOCIATION 

For  every  track  pairs  within  an  association  gate,  let  N  denote  the  total  number  of  common  track  points 
and  let  M  denote  the  total  number  of  times  both  the  tracks  were  updated  simultaneously,  so  that  N-M 
equals  total  number  of  times  at  least  one  track  was  not  updated.  Then  the  test  statistic  tor  evaluating 
how  close  the  two  tracks  are,  over  the  entire  track  matching  interval,  can  be  obtained  by  constructing 
the  closeness  score,  CSCt^)  defined  as: 
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M  .  1  -  •  -  1 

CS(tM)=X(Xkt|c/t„)-X'{tk/tk))’'(P'(tk/ti)  +  pJ(H/tk)r'(X'(tk/tt)-X'(H/tk))(«) 

k=l 

which  is  a  cumulative  sum  over  M  track  match  points.  The  probability  distribution  of  correct 

association  of  two  tracks  X*  and  X*^ ,  denoted  by  Pca  is  then  obtained  by  assuming  independence  of 
the  error  between  the  two  track  estimates  and  by  testing  against  the  threshold 

corresponding  to  the  chi-square  distribution  with  nM  degrees  of  freedom: 

P„  =Pr(CS(tM)S^aM/H„)  =  l-a  H-4) 

In  order  to  incorporate  the  dependence  between  the  track  estimates  because  of  the  common  process 
noise,  the  closeness  score  is  modified  as: 


cs'(tv,)  =  l!xntk 

k=l 


/tk)-Xntk/tk)}^P^(tk/tk)‘'(X‘(tk/tk)-Xj(tk/tk)} 


(4-5) 


Let  CS,(tM)  *  0  CS*(tM)/  N,  where,  0  is  a  sensor  dependent  constant.  Then  at  every  time  update,  the 
closeness  score  is  converted  to  a  Figure  of  Merit  (FOM)  by  mapping  it  to  a  range  [0,1]  by  means  of 
the  following  transformation; 


FOM(tM) 


1 

l  +  CSjCt^) 


b 

b  +  CS*(tM) 


(4-6) 


where,  b  =  N/  0.  This  FOM  is  also  used  as  a  goodness-of-fit  statistic  for  track-to-track  association. 
Incorporating  the  cross-covariance  in  the  modified  closeness  score,  the  probability  distribution  of 
correct  association  denoted  by  P*ca  becomes 


Pca  ~  Pr(CS  (t^ )  ^  ^cs‘,nM  ^  ^o)  “  ^ 


(4-7) 


II 


A  measure  of  effectiveness  which  was  found  useful  for  this  track-to-track  association  problem  is  the 
FOM,  given  in  (4-6).  Corresponding  to  this  FOM,  the  probability  of  correct  association^  P*ca 
cumulative  probability  distribution  function  given  by: 


PcA  ®  Pr(FOM  ^  XFOM.nM) 

b 


=  Pr( 


b  +  CS  (tv,) 


<X 


FOM.nM 


) 


(4-8) 


where,  f^j-  ^  (z)  is  the  chi-square  density  with  nM  degrees  of  freedom,  and  nvi  given  by 


'cs\nM 


=  b(: 


-FOM.nM 


-1) 


(4-9) 


It  is  to  be  noted  that  numerical  computation  for  P*ca  nsing  (4-8)  is  complicated.  Instead,  an  indirect 
approach  using  (4-9)  is  preferable  because  it  lends  itself  to  the  use  of  readily  available  tables  for  chi- 
square  densities.  Using  such  a  table,  trade-off  studies  for  a  goodness-of-fit  test  can  be  readily 
performed,  as  shown  in  table  1. 
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Table  1.  Probability  of  Correct  Association  Trade-Off  Study 


PCA 

M 

^S*.nM 

b 

^FOM.nM 

2 

6.30 

119.70 

0.95 

623.70 

0.99 

Pr(FOM<A.FOM.aM)  =  0.90 

J 

10.9 

207.10 

0.95 

1079.10 

0.99 

5 

20.60 

391.40 

0.95 

1 

2039.40 

0.99 

2 

5.23 

99.37 

0.95 

517.70 

0.99 

Pr(FOM  <  XpoM.nM)  ~  0-95 

3 

9.40 

178.60 

0.95 

930.60 

0.99 

5 

18.50 

351.50 

0.95 

1831.50 

0.99 
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This  table  is  useful  for  a  system  designer  who  may  choose  to  select  the  level  of  Pca  a-priori.  For 
example,  if  it  is  desired  to  have  Pqa  ~  0-95,  then  (4-8)  is  used  to  fix  the  threshold,  .  for  each 

M,  which  corresponds  to  a  level  of  significance  of  0.05.  In  this  study,  it  is  assumed  that  the  three 
sensors  are  equipped  with  six  dimensional  Kalman  filters.  Hence,  n  =  6.  At  the  end  of  every  fast 
track  update  time  t^,  the  slower  tracks  are  predicted  to  t^.  The  predicted  tracks  are  subjected  to  a  chi- 
square  test  for  association  using  si.x  degrees  of  freedom.  Every  time  the  predicted  track  passes  this 
test,  the  counter  for  N  is  incremented  by  1  and  the  number  of  track  match  points  M  is  adjusted.  .\t 
time  t^  if  the  number  of  track  matching  points  increase,  then  the  counter  for  M  is  incremented  by  1, 
otherwise,  it  is  left  unchanged.  Using  this  value  of  M,  for  each  fast  track  i  and  every  slow  track  j,  the 
closeness  score,  CS^Ct^),  is  computed  using  (4-5).  For  each  t^,  these  values  of  CS*(t|v))  are  compared 
with  which  can  be  obtained  from  a  look-up  table.  Final  selection  for  association  is  made  tor 

the  track  pair  which  has  the  highest  level  of  CS*(tM)  and  exceeds  • 

Table  I  also  shows  a  trade-off  between  b  and  ^FOM,nM*  addition,  the  choice  of  b  affects  the 
parameter  0  as  shown  in  (4-6).  But  the  choice  of  9  is  dependent  upon  the  sensor  mode  and  other 
parameters.  Hence,  there  is  a  trade-off  between  ^FOM,nM  9  as  shown  in  table  I . 

4.2  PROBABILITY  OF  FALSE  CORRELATION 

Selection  of  the  threshold  for  the  test  of  association  has  direct  bearing  upon  the  probability  of  false 
correlation  Ppc-  Analytical  derivation  of  this  expression  in  closed-form  follows  closely  the 
methodology  developed  in  [8].  In  this  analysis,  it  is  assumed  that:  a)  false  tracks  are  uniformly 
distributed  over  the  six  degree  of  freedom  chi-square  gate  for  association,  and  b)  false  correlations  are 
independent  tor  each  update  which,  essentially,  implies  random  clutter. 

Let  denote  the  size  of  the  six  degree  of  freedom  chi-square  gate  for  association.  This  gate  size 
must  be  selected  to  let  sufficient  number  of  slow  tracks  through  this  window.  Assuming  that  a  false 
track  is  uniformly  distributed  within  this  gate,  the  cumulative  probability  distribution  function  of 
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y(tk)=  [x‘(tk/tk)-xj(tk/tk)f(P^(tk))’'[x‘(tk/tk)-xj(tk/tk)]  (4-10) 

is  given  by 


Pr[y(WSX]=(:^)2  forOSXSX, 


(4-11) 


thus,  the  pdf  of  y(tk)  is 


f[y(tk)]  = 


n(y)^ 


for  0  <  A.  <  Xg 


(4-12) 


Assuming  independence  of  false  association  at  each  the  joint  pdf  of  M  false  correlations  is  the 
convolution  of  M  pdfs  given  by  (4-12).  Hence,  given  a  test  threshold  Xpcj  the  false  track  correlation 
probability  is 


Ppc  =  Jg(x)dx  (4-13) 

0 

where,  g(x)  is  the  convolution  of  M  pdfs  given  in  (4-12).  Evaluation  of  (4-13)  in  closed-form  is 
analytically  intractable  because  of  the  nature  of  f(y(t|t))  given  in  (4-12).  For  this  reason,  g(x)  is 
approximated  by  a  Gaussian  pdf  whose  mean  and  variance  are  given  below 


^  =  ME[y]  = 


MnA.g 
n  +  2 


and 


1 

(7 


=  MVar[y]  = 


4MnA.^ 

(n  +  4)(n  +  2)^ 


(4-14) 


(4-15) 
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Therefore,  the  false  track  correlation  probability  is  approximated  by 


-00 


TIT 


where, 


Xpc  - 


^cs’.nM  ” 


(4-16) 


(4-17) 


Based  on  this  approximation,  table  2  presents  results  of  trade-off  study  for  n  =  6,  which  involves  Pca> 
X„ ,  M,  and  Ppc-  This  table  shows  that  for  a  chosen  Xg,  Ppc  decreases  as  M  (the  number  of  track 
matching  points)  increases,  for  a  desired  level  of  Pca-  On  other  hand,  if  it  is  desired  to  hold  P pc  to 
a  low  value,  then  either  M  must  be  increased,  which  will  delay  decision  on  track  pair  association  and 
may  cause  track  staleness,  or  Pqa  must  be  decreased.  The  same  results  could  also  be  achieved  by 
making  the  gate  window  Xg  tighter  which  would ,  in  turn,  reject  good  candidate  tracks  for  association, 
and  force  the  system  designer  to  reduce  Pca- 

It  is  to  be  noted  that  the  analysis  presented  in  this  section  is  an  analytical  tool  that  can  be  used  ott- 
line,  to  predict  the  Measure  of  Effectiveness  (MOE)  of  a  track  fusion  algorithm  without  resorting  to 
elaborate  Monte  Carlo  simulation.  For  this  reason,  the  algorithms  for  computing  probabilities  of 
correct  association  and  felse  correlation  have  not  been  implemented  in  the  software  developed  in  this 
contract. 
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Table  2.  Probability  of  False  Correlation  Trade-Off  Study 


PCA  = 

0.90 

PCA  = 

'0.95 

M 

X 

o 

s 

X 

g 

90% 

95% 

90% 

95% 

Ppc 

Ppc 

Ppc 

Ppc 

I 

0.90 

0.68 

0.99 

0.90 

2 

0.81 

0.46 

0.96 

0.73 

0.72 

0.29 

0.92 

0.54 

5 

0.54 

0.10 

0.81 

0.27 
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SECTION  5 


KINEMATIC  FUSION  OF  ASYNCHRONOUS  TRACKS 

In  this  section,  the  algorithm  for  Kinematic  Track  Fusion  (KTF)  is  described.  The  KTF 
algorithm  is  basically  a  weighted  covariance  fusion  algorithm  where  the  candidate  tracks  for 
fusion  are  weighted  by  certain  covariances.  The  weighted  tracks  are  then  summed  to  obtain  the 
fused  track  as  shown  in  (2-9)  -  (2-10).  It  is  assumed  that  targets  are  detected  by  three  sensors 
(Radar,  IR,and  Laser  Radar).  Details  of  the  sensor  models  are  described  in  Appendix  A.  After 
detection,  the  sensors  perform  report  to  track  association  and  transmit  these  tracks  to  a  central 
station  by  communication  links  (e.g.,  TADIL  A,  B,  or  J).  It  is  assumed  that  the  Radar  transmits 
the  tracks  at  a  rate  of  4  sec.,  the  IR  sensor  transmits  at  a  rate  of  2  sec.  and  the  Laser  radar 
transmits  at  a  rate  of  1  sec.  These  sensors  have  non-unity  probability  of  detection  (P □<  1). 
Hence,  the  tracks  created  after  performing  report  to  track  association,  may  contain  gaps  due  to 
missed  detection.  In  addition,  false  tracks  may  be  created  which  will  result  in  the  number  of 
tracks  being  greater  than  the  number  of  targets.  These  tracks  are  processed  at  the  central  station 
where  track  synchronization  (described  in  Section  3)  and  association  (described  in  Section  4)  are 
performed.  Tracks  are  associated  by  minimizing  closeness  score,  and  the  track  pair  that  gives 
rise  to  minimum  closeness  score  are  then  fused  using  the  Kinematic  Track  Fusion  (KTF) 
algorithm  described  in  this  section.  Track  fusion  is  useful  for  airborne  surveillance  and  C^I, 
because  two  sensors  tracking  the  same  target  may  contain  different  kinds  of  information  about 
the  target 

It  is  well  known  [2]  that  when  two  sensors  track  the  same  target,  the  tracks  become  correlated 
due  to  target  maneuver  noise  which  is  the  same  for  the  two  trackers  tracking  the  target.  It  has 
also  been  shown  that,  incorporation  of  cross-correlation  between  tracks  in  the  track  fusion 
algorithm  (3]  results  in  reduction  of  variance  of  position  and  speed  which  are  assumed  to  be  the 
states  of  the  two  track  filters.  However,  published  results  in  open  literature  consider  fusion  of 
synchronous  tracks  that  are  created  by  perfect  sensors  of  Pq  *  1.  The  KTF  algorithm  described 
this  section,  consider  fusion  of  asynchronous  tracks  that  contain  missing  data  (introduced  by 
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sensors  that  have  non-unity  probability  of  detection).  An  efficient  algorithm  for  computing 
cross-correlation  between  these  tracks  is  derived.  At  any  time  t^,  if  a  track  from  one  sensor  is 
not  updated,  then  this  track  is  propagated  to  t|(  from  its  last  update  at  tjj.,  and  cross-correlation 
between  this  propagated  track  and  tracks  from  other  sensors  are  computed.  For  asynchronous 
tracks,  this  cross-correlation  is  computed  by  modifying  the  equation  for  cross-covariance  given 
in  (2-8).  Then  the  fused  track  and  its  covariance  is  computed  by  incorporating  the  cross¬ 
covariance  in  the  track  fusion  algorithm  as  shown  in  (2-9)  and  (2-10). 

5.1  DESCRIPTION  OF  TARGET  AND  TRACK  MODELS 

In  this  section,  it  is  assumed  that  four  targets  are  moving  at  the  same  speed  of  0.3 1  km/sec  at  a 
constant  altitude  of  9. 10  km.  The  target  dynamic  is  modeled  by  means  of  a  state  variable  model 
which  has  six  states  (three  dimensional  positions  and  speed).  Target  maneuver  noise  is  modeled 
by  a  three  dimensional  white  Gaussian  noise  vector.  This  target  model  is  updated  every  one 
second  or,  AT  =  1 .  Sensor  accuracy  is  modeled  by  white  Gaussian  noise  vectors.  All  the 
sensors  are  assumed  to  be  colocated  at  the  same  altitude  as  the  targets. 

At  the  fusion  center,  tracks  from  one  sensor  are  associated  with  tracks  from  the  other  sensor,  one 
at  a  time.  Let  the  tracks  from  one  sensor  be  denoted  by  i  and  tracks  from  the  other  sensor  be 
denoted  by  j.  Tracks  from  sensor  i  (or  j)  are  classified  as  follows: 

Updated  tracks  from  sensor  i  at  time  t^+i  are  denoted  by  i,  and  are  counted  as  NIU(t|j+,).  These 
updated  tracks  are  further  subdivided  as  follows: 

•  ij  tracks  updated  at  ti^  and  at  ti^+i  are  counted  as  NlUlft^+i)  and  labeled  as  I 

•  i|  tracks  not  updated  at  t|(  but  updated  at  ti^^i  are  counted  as  NIU2(t|t+|)  and  labeled  as  2 

•  newly  acquired  updated  tracks  at  t|j+|  are  counted  as  NIU3(t|(+i)  and  labeled  as  3 

After  a  new  track  is  received  at  the  fusion  center,  it  is  correlated  with  tracks  from  the  other 
sensor  and  then  added  to  the  existing  NIUl(tk,.i)  track  file.  Separate  processing  of  newly 
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acquired  track  is  necessary  because  the  initial  value  for  the  recursive  cross-covariance  matrix  for 
this  track  is  zero.  Once  the  cross-covariance  matrix  is  initiated,  the  newly  acquired  track  can  be 
classified  as  an  i|  track.  Hence, 

NIU(tk)  =  NlUl(tk)  +  NIU2(tk)  +  NIU3(tk) 

Non-updated  tracks  at  t^+i  are  denoted  by  i2  and  are  counted  as  NIXU(t|{>,).  These  non-updated 
tracks  are  further  subdivided  as  follows: 

•  i2  tracks  updated  at  t^  but  not  updated  at  ti^*]  are  counted  as  NIXU I  (t|(+ 1 ) 

•  i2  tracks  not  updated  at  t^  and  not  updated  at  t^+i  are  counted  as  NIXU2(tj5,.i) 

Non-updated  NIXU2  tracks  are  deleted  if  they  are  not  updated  for  ten  successive  update  times. 
Otherwise,  computer  resources  necessary  to  implement  this  algorithm  will  be  excessive. 
Furthermore,  because  of  the  recursive  nature  of  the  fusion  algorithm,  it  is  necessary  to  store 
track  statistics  only  for  two  successive  time  updates  which  reduces  the  burden  of  computation. 
Hence,  total  numbers  of  non-updated  tracks  at  any  time  ti^.^!  are: 

NIXU(Vt)  “  NIXUl(tk^,)  +  N1XU2(V,) 

Format  of  a  track  file  from  sensor  i  is  given  below: 

i(track  #,  track  id(lor2or  3)  at  current  timeftk,.,),  track  id(lor2)  at  t^,  y 

Since  the  tracks  are  asynchronous,  let  the  update  rate  for  tracks  from  sensor  i  be  denoted  by  NI 
and  the  update  rate  for  tracks  from  sensor  j  be  denoted  by  NJ.  Let  the  minimum  update  rate  be 
denoted  by  NMIN  -  Min(NI  W-  Tracks  from  these  sensors  are  used  to  perform  Kinematic 
Track  Association  (KTA)  by  minimizing  closeness  score,  and  the  tracks  which  give  rise  to 
minimum  closeness  score  are  fused  by  means  of  a  Kinematic  Track  Fusion  (KTF)  algorithm. 


5.2  DESCRIPTION  OF  WEIGHTED  COVARIANCE  KTF  ALGORITHM 


At  time  t^,  tracks  from  the  i  and  j  sensors  may  have  different  classification  according  to  the  state 
of  update  (Track  ID  -  1  if  it  is  updated  at  t^  and  at  t^.i.  Track  ID  >  2  if  it  was-updated  at  t|(.|  but 
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is  not  updated  at  ti^  and  Track  ID  =  3  if  it  is  a  newly  acquired  track).  Hence,  the  algorithm  for 
KTA  must  incorporate  this  information  at  every  update  time.  Since  the  KTA  algorithm 
described  in  Section  4,  uses  cross-covariance  information  for  track  matching,  cross-correlation 
between  tracks  of  different  classification  are  computed  separately.  Derivations  of  these  cross- 
correlation  matrices  are  described  below.  Let 

NI 

QI  =  J;<|)P-'(1)G(1)QGT(1)((|)P-'(1))^  (5-1) 

P=l 

NJ 

QJ  =  2;'I>‘’''(1)G(1)QG^(1)((|)P-'(1))^  (5-2) 


Q2=  2I<l»‘’‘'G)G(l)QG^(l)((i)P-'(l))^ 


and  let  the  recursive  equation  for  cross-covariance  be  initialized  at  time  to  as 

(tfl  ’  ^0 i ^0  >  ^0 )  =  0  •  For  a  newly  acquired  track  at  time  t^,  the  initial  cross-covariance  at  t^.j 
will  also  be  zero. 

Let  the  tracks  from  sensor  i  be  the  faster  tracks  which  are  indexed  by  N3  and  let  the  slower 
tracks  from  sensor)  be  indexed  by  N,.  Then  - 1^  =  AT  (N,)  NI  =  AT  (NJI)  NJ  where, 

N3  =l,2,....,NJI.  To  illustrate,  the  tracks  from  Laser  Radar  are  updated  every  1  sec.,  or  NJ  =  1, 
and  the  Radar  tracks  are  updated  every  4  sec.,  or  NI  =  4.  Hence,  NJI  =  4  and  tracks  from  both 
the  sensors  are  updated  simultaneously  every  4  sec.  For  N3  =l,2,....,NJI  -  1,  the  Radar  tracks  are 
predicted  every  NJ  =»  NMIN=1  sec.  as  follows; 


F'''(^(N1-1)+N3  ^  ^Nl-|)  “*I*'^'^^F^*(^(NI-I)+(N3-I)  ^  ^Nl-1)  ■*’Q2 

It  is  to  be  noted  that  for  N3  =  NJI, 
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^**(^(Nl-l)+N3  ^  *(NM)+N3)  “  ^(^(NM)+N3)  ®  P  ~  ^  P(NI-I)+N3)^  P(NI-l)+N3)l*l* 

r^'*P(NI-l)+N3-I  ^  *(Nl-l)+N3-l)  “  ^P(Nl-l)+N3-|)]  ^  P(>4I-I)+N3)^  P(NI-1)+N3) 

NI 

-  P~K'*(t(^^|.|)^N3)^'*P(NI-l)+N3)]{S^**  GW(t(|^l.l)+N3_p)} 

p-l 

(5-6) 

and, 

^^'(^(Nl-l)+N3  ^  ^NI-1)h-N3)~^(*(NI-I)+N3)  =  P  “  ^^*(^('NI-I)+N3)^^  (^(NI-I)+N3)1‘{> 

r^^*(^(Nl-l)-t-N3-l  /^■NI-I)+N3-l)”^(^CNI-l)+N3-l)l'^  (t(Nl-I)+N3)^  (t(NI-l)+N3) 

NMIN 

P-) 

(5-7) 

Since,  t(Ty„.i)+NJi  ®  *Ni.  the  cross-covariance  matrix  between  X''(tNi  /t^i)  and  X^'(tNi  /tNi)is 
computed  as, 

P‘'’^'(tNl.tNi;tNl.tNl)=P  -K‘'(tN,)H“(tN,)]<l>^'P“-^'(tNl.NJ>tN,.Nj;t  Nl-NJ«t(Nl-l)) 

(^,NMIN)T  p.Kjl(,  )Hf(tN|)r+[I  -K"(tN,)H''(tN,)lQ2tI-K''(tNl)Hj'(tN|)r 

(5-8) 

This  recursive  equation  is  initialized  by  P‘''^'(tNi-NJ  ’^ni-nj  ■'t!sji_Nj>tNi-i)  which  was 
computed  at  tN,.Nj.  Now,  using  (2-9)  and  (2-10)  the  fused  track  is  computed  as, 

P^(tNl  /  t^l)  *  P''(tNl  ^  tNl)''’P^*(tNl  ^  ^Nl)”  P''’'''(tNl>tNlitNl»tNl)  ~  P'*’^'(tNl»tNi;tNl>tNl) 

X‘'(tNl/tNl)*X“(tN,/tN,)  +  (P“(tN,/tN,)-P“‘j'(tN„tN,;tN„tN,)](P^)‘'‘. 

[Xj‘(tN,/tN,)-X‘'(tN,/tN,)l 

P^t^,  /tN,)»P‘‘(tN,/tN,)-[P“(tN,  /tN,)-P“’^'(t^j,,tN,;tN,.tN,)l(P^)-'. 

[P"(tN!  /tNl)-P‘'’^'<tNl.tNi;tNl.tNl)j’' 

(5-9) 

Similarly,  when  the  fast  track  i,  has  been  updated  and  the  slow  track  j  i  has  been  predicted  to  the 
update  time  of  the  fast  track,  cross-covariance  between  X'*  (t(Ni-i)+N3  ^  t(Ni-o+N3) 
Xj'(t(N,.,)^N3/tN,.,)  forNj=  1,2 . NJl-l  can  be  derived  as. 
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P'*’‘'*(^NM+N3»^NI-l+N3’^Nl-l+N3’^Nl-l)  =  (tNl-l+N3)^'* (^Nl-l+N3)]‘l> 

P'*’^*(t(Nl-l)+N3-l  >  ^  (N 1-1  )+N3-l  i  *^(N  1-1  )+N3-l  ’  t  (N  I-l )  X*!^ 

[I  -K"(tN,.KN5)H"(tNI.|-N))]Q2' 


NMIN^T 


(5-10) 


Now  at  every  N3,  the  fast  track  update  index,  the  fused  track  and  its  covariance  can  be  obtained 
by  using  (2-9)  and  (2- 1 0)  similar  to  (5-9). 


Cross-correlation  between  the  ij  tracks  X''(t(Ni-i)+N3  ^  ^ni-I)+N3)  which  are  updated  at 

t(Ni-i)+N3  and  were  last  updated  at  and  the  j,  tracks  Xj'(t^,_,)+tM3  ( ^Ni-i))  <^an  be 

derived  as. 


(^Nt-l+N3’tNl-l+N3itNl-l+N3*^l-l)  -  [I  '  X 


'(t 


Nt-t+N3 


)H''(tNl-l+N3)M>'^*- 


P'^’'''(^NI-1)+N3-1  ’^(Nt-t)+N3-Mlit(Nl-t)+N3-l’t(Nl-l)X<l> 

[I  -  K‘'(tj^l.H.N3)H''(tNl-l+N3)]Q2 


NMIN)T^  (5-11) 


This  equation  is  initialized  by  P'^^‘(t(Nt-i)+N3-i»t(Nt-i)+N3-Miit(Ni-i)+N3-i  -tNi-i)  '^^ich  was 
computed  at  t(N,.,)+N3,,.  Computation  of  fused  tracks  is  now  performed  using  (2-9)  and  (2-10). 
Cross-correlation  between  i;  tracks  X'^(t(^i_i)+)v43  /  t(Tsj|_|)^.N3.i)  counted  as  NIXUl  which 
were  last  updated  at  t(N|.i)+N3.i  and  ji  tracks  X^*(tQvjj_|)+jsj3  /  t(|sji_i))  is  described  below: 


-t+N3 » t  N 1-1+N3-1  ’  ^  N  t.l+N3 »  ^Nl-l )  =  P'*’’’ '  (*(Nt-l)-*-N3-t  >^(Nl-t)+N3-li^(Nl-t)+N3-l» 

t(Nl-t)X<l»''^"^)^ +Q2 


(5-12) 


This  equation  is  initialized  by  P'‘’-''(t(Ni-i)+N3-i.t(Ni-n+N3-i;t(Ni-i)+N3-i-tNi-i)  which  was 
computed  at  t(Ni.i)*N3-i-  Computation  of  fused  tracks  is  now  performed  using  (2-9)  and  (2-10). 
Cross-correlation  between  ii  tracks  X'^(t(^i_i)^.N3  /  t(Ni_i).,.;vj3.[^2)  counted  as  NIXU2  which 
were  last  updated  at  t(Ts(|.,)+N3.M2  and  j,  tracks  X'’*(t(N|_|)+{s(3  ^  is  described  below; 
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P'^'**(^Nl-l+N3»*Nl-t+N3-M2’^NI-l+N3>^NM)  P  '■*  (t(NI-l)+N3-l»^(NI-l)+N3-M2’^(NI-l)+N3-l’ 

»  VA  NMfN  \T  . 

(5-13) 

This  equation  is  initialized  by  P'^^'  (t(Ni-i)+N3-i » t(Ni-i)+N3-M2  ♦  t(Ni-i)+N3-i » ^Ni-i )  which  was 
computed  at  t(Ni.i)+N3.i.  Computation  of  the  cross-correlation  between  tracks  and  i,  as  well  as 
U  tracks  is  performed  using  identical  methodology  and  further  details  are  omitted.  It  is  to  be 

noted  that  the  j,  tracks  XJ^t(^,_,)^N3  /  t,N,_2))  counted  as  NJXUl  were  j,  tracks  at  t,N,.2). 
These  tracks  are  propagated  to  t(N|.i)+N3  and  then  correlated  with  i2  tracks.  The  i2  tracks  are 
processed  similarly.  To  illustrate,  cross-correlation  between  ii  tracks  counted  as  NIXUl  and  j2 
tracks  counted  as  NJXUl  is  given  by, 


P'^^*(*NI-I+N3’^NI-1+N3-iUnI-1-^N3’^I-2)  ®  *1^  *P  (tcNi.|)+N3.i .t(^l.l)+N3-l ’ ^fNl-l)+N3-! ’ 


^(NI-2)Xt  ) 


(5-14) 


This  equation  is  initialized  by  P'*’-’*(t(>ji.|)+N3.i>*(NM)+N3-ii^fNi-i)+N3-i*^Ni-2)  which  was 
computed  at  t(Ni.i)+N3-i-  Similarly,  cross-correlation  between  ia  tracks  counted  as  NIXU2  which 
were  last  updated  at  t(Ni.,)4.N3.M2  and  j:  tracks  counted  as  NJXUl  is  given  by, 


P‘^''^(^Nl.|)+N3**(NW)+N3-M2it(NH)+N3**Nl-2)  *'^^*P'^'’'(*(NH)+N3.|»*(N!-l)+N3-M2i 

^(NI-l)+N3-l**(NI-2)X+^*'^)^  +Q2 

(5-15) 

This  equation  is  initialized  by  P‘^^'(t(Ni_,)+N3-i**(Ni-i)+N3-M2i^(Ni-i)+N3-i*^Ni-2)  which  was 
computed  at  t,N,.|,.,.N3.i.  Cross-correlation  between  ii  tracks  counted  as  NIXUl  which  were  last 
updated  at  t,Ni.,,*N3.i  and  j2  tracks  counted  as  NJXU2  is  given  by. 


pi2.j2^ 


^NM+N3*^(Nt-l)+N3-l>^Nl-l+N3»^l-M3)  “ ^  ’^^(t(Nl-l)>N3-l •*(Nl-l)+N3-l > ^Nl-I)+N3-I 

tfN»-M33  X^  +  Q2 


(5-16) 
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This  equation  is  initialized  by  P'‘’-'^(t(Ni-i)+N3-i»t(Ni-i)+N3-i;t(Ni-i)+N3-i’tNi-M3)  which  was 
computed  at  t(Ni-i)+N3-i-  Similarly,  cross-correlation  between  i2  tracks  counted  as  NIXU2  which 
were  last  updated  at  t(ig|.,)+N3.M2  and  j2  tracks  counted  as  NJXU2  is  given  by, 


^(Nl-M3))('t’ 


NMIN)T^q2 


(5-17) 


For  N3  =  NJI,  slower  tracks  from  sensor)  are  updated  and  cross-correlation  between  these 
updated  tracks  and  the  updated  tracks  from  sensor  i  is  computed.  For  i]  tracks  counted  as  NIUl 
and),  tracks  counted  as  NJUl,  the  cross-correlation  matrix  is  shown  in  (5-8).  Cross-correlation 
between  tracks  of  other  classifications  (i.e.,  i2  and  j,,  i2  and  j2,  and  i2  and  j2)  can  be  derived  using 
the  methodology  described  earlier. 


Initialization:  Three  track  updates  from  each  sensor  are  required  to  initialize  the  KTF  algorithm. 
Until  the  third  track  update  from  the  slower  sensor  is  received  at  the  fusion  center,  tracks  from 
the  faster  sensor  are  monitored  and  the  track  labels  are  stored.  Labels  of  i|  tracks  that  have  been 
updated  at  the  current  time  update  t)j+,  and  the  previous  time  update  t^  are  stored  as  1 .  The  i| 
tracks  that  have  been  updated  at  the  current  time  update  t^*,  but  were  not  updated  at  are 
labeled  as  2.  These  i,  tracks  were  created  earlier  but  contain  missing  data.  The  i2  tracks  that 
were  created  earlier  but  have  not  been  updated  recently  are  predicted  one  step  ahead. 

This  algorithm  was  coded  in  MATLAB  and  the  simulation  results  are  shown  in  Section  6.  An 
alternative  algorithm  based  on  information  fusion  was  also  implemented.  Details  of  this 
software  are  described  in  the  next  section. 


5.3  DESCRIPTION  OF  INFORMATION  FUSION  APPROACH  TO  KTF 

Recently,  there  have  been  considerable  interest  in  the  information  fusion  approach  [9-11]  to 
centralized  KTF.  This  approach  is  different  from  the  weighted  covariance  fusion  described  in 
Section  2.  The  difference  is  in  the  methodology  employed  to  incorporate  the  track  correlation 
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information  contained  in  the  remote  tracks  that  are  transmitted  to  the  fusion  center.  In  the 
weighted  covariance  track  fusion  methodology,  the  algorithm  assumes  that  the  remote  tracks  are 
correlated  over  time  and  are  also  correlated  because  of  the  input  maneuver  noise  which  is 
common  to  both  the  tracks.  The  KTF  algorithm  described  in  Section  2,  computes  cross¬ 
covariance  between  the  tracks  recursively  and  updates  it  as  soon  as  new  information  is  received 
at  the  fusion  center.  Then  this  algorithm  weights  the  tracks  by  incorporating  the  cross¬ 
covariance  matrix  and  the  weighted  tracks  are  fused  using  (2-9).  On  the  other  hand,  the 
information  fusion  approach  recognizes  the  fact  that  remote  tracks  are  correlated  in  time  and  the 
measurement  noise  could  be  correlated.  But  this  algorithm  does  not  take  into  account  cross- 
correlation  between  tracks. 


In  this  section,  it  is  shown  that  the  KTF  algorithm  (2-9)-(2-10)  which  was  derived  in  [3]  is  not 
exact  but  an  approximation  of  the  optimal  solution.  In  feet,  a  published  study  [15]  showed  that 
the  performance  of  the  fusion  algorithm  derived  in  [3]  (referred  to  as  the  state  vector  fusion 
method”)  is  consistently  worse  (biased)  than  the  optimal  method  (referred  to  as  the 
‘'measurement  fusion  method”).  It  was  pointed  out  that  the  method  of  combining  tracks  using 
state  vector  fusion  is,  in  general,  suboptimal.  However,  this  study  failed  to  recognize  that  the 
real  reason  behind  the  performance  degradation  of  this  algorithm  is  an  error  in  its  derivation.  In 
this  section,  the  nature  of  the  approximation  made  in  [3]  is  examined.  Performance  of  this 
approximate  algorithm  is  compared  with  that  of  a  track  fusion  algorithm  based  on  information 
fusion  in  the  next  section. 

The  KTF  algorithm  (2-9)-(2-10)  derived  in  [3]  is  based  on  the  following  equality 

E[x^|D‘]  =  x’  (5-18) 

where.  D'  denotes  the  information  (accumulated  measurements)  from  sensor  i.  This  equality  is 
incorrect  because 
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E[xj|D‘]=  Jxjp(Dj|D‘)dDj  =  |E[x|Dj]p(Dj|D‘)dDj 
=  J[[xp(x|Dj)dx]3(Dj|D‘)dDj 
=  Jx|p(x|Dj)p(Dj|D‘)dDjdx 
;6  JxJp(x|D‘,Dj)p(DjlD‘)dDjdx 
=  Jxp(x|D‘)dx 
=  E[x|D‘]  =  x‘ 

Hence,  it  is  clear  that  (5-18)  will  hold  if  and  only  if  p(x|D-’)  =  p(x|D',D-')  which  is  incorrect. 
Therefore,  (5-18)  can  only  be  considered  an  approximation  and  so  is  the  KTF  algorithm  (2-9)- 
(2-10). 


An  alternative  KTF  algorithm  based  on  the  concept  of  information  fusion  [9]  is  described 
below; 


x‘'(t 

P''(t 


k^,/t„^,)  =  P''(tk.,/tk^,)(P‘(tk^,/tk^,)-'x'(tk^,/tk^,)  +  pj(tk^,/tk^,)-‘xj(tk^,/tk.,) 
-P'(tk+i/tk)'‘x‘‘(tk.i/tk)-pj(tk^,/t„)-'xj(tk^,/tk)  +  P^tk^,/tk)-'x‘'(tk.,/tk) 

k.,  /tk.l)"'  =  p‘(tk-l  /tk.,)’'  +p'(tkH  /tk.,)-'  -Pktk^l  /tk)"'  -p'ok.l  /tk)“'  + 
p"(tk.l/tk)-‘ 


(5-20) 


This  algorithm  is  exact  if  the  two  sensors  communicate  each  time  (full-rate  communication) 
they  receive  measurements  and  i^Kiate  tracks  or  when  the  process  noise  is  zero  (deterministic 
system).  Hence,  in  general,  this  algorithm  is  also  approximate.  It  is  to  be  noted  that  the  time 
correlation  information  contained  in  the  fusion  algorithm  (5-20)  is  removed  by  subtracting  out 
the  prior  information  contained  in  the  last  terms  of  this  equation.  In  addition,  this  algorithm 
does  not  take  into  account  the  cross-correlation  due  to  input  maneuver  noise.  In  order  to 
compare  the  performance  of  the  two  algorithms,  the  fused  covariance  equation  in  (5-20)  is 
rewritten  as. 
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(5-21) 


i-i 


P‘^=(P'  +PJ  -P'  -PJ  +P''  ) 


=  P' 


=•:  1.  =;F“* 


■^r 


p'+pj_pj(p'  +pj  -p^  )P' I  pj 

= p‘  -  p‘(p' + pj  -  p*)''(p‘  -  p'') 

and, 

P#  =  pj(pi"'  +pr'  -pp"’)p‘ 


where,  P' ,  pj ,  and  P^  denote  the  predicted  covariances,  and  for  the  sake  of  simplicity,  time 
dependence  has  been  omitted.  It  can  be  seen  from  (5-2 1 ),  that  one  reason  why  the  weighted 
covariance  fusion  algorithm  (2-9)-(2- 1 0)  is  not  exact  is  due  to  the  fact  that  no  prior  information 
is  used  in  this  formulation.  It  should  also  be  noted  that  when  the  tracks  from  the  two  sensors  are 

assumed  to  be  independent,  then  the  cross-correlation  is  zero,  or,  P*^  =  P**  =  0 .  and  the  two 
fusion  algorithms  become  identical.  Furthermore,  when  the  two  sensors  transmit  tracks 
asynchronously,  (5-20)  also  becomes  approximate.  This  is  so,  because  the  assumption  that  the 
two  sensor  measurement  errors  are  independent,  becomes  invalid  due  to  the  propagation  of  the 
process  noise  through  predicted  covariances[9]. 


This  algorithm  was  coded  in  MATLAB  and  details  of  this  software  are  described  in  the  next 
section. 


28 


SECTION  6 


SIMULATION  RESULTS 

This  section  describes  the  results  of  simulation  of  the  two  algorithms  described  in  Section  5.  A 
three  sensor,  four  target  scenario  shown  in  Figure  1.  was  simulated  using  MATLAB  4.0 
software  package  and  a  supporting  platform.  In  this  scenario,  three  sensors  (Radar,  IR.  and 
Laser  Radar)  are  colocated  and  mounted  on  a  surveillance  platform  which  is  stationary  at  an 
altitude  of  9. 10  km.  Threat  aircraft  altitude  is  also  at  9. 10  km.  Speed  of  the  target  aircrafts  are 
all  the  same  at  0.3 1  km/sec.  Threats  are  separated  in  the  y-axis  by  37. 1 0  km.  Separation 
between  sensor  platform  and  threats  is  371  km.  Details  of  the  sensor  models  are  given  in 
Appendix  A.  For  this  simulation,  a  sensor's  pointing  is  with  respect  to  North  and  azimuth  is 
measured  clockwise.  Monte  Carlo  simulation  is  performed  to  generate  tracks  from  each  sensor. 
These  tracks  are  used  at  the  fusion  center  by  the  two  algorithms  (covariance  weighted  fusion  and 
information  fusion)  to  perform  track  to  track  fusion. 

The  weighted  covariance  fusion  algorithm  can  be  executed  through  a  Graphical  User  Interface 
(GUI)  which  is  implemented  in  Visual  Basic.  Parameters  of  each  sensor  can  be  chosen  by  a  user 
by  entering  the  values  in  the  screen  for  sensors,  or  by  accepting  the  default  parameters.  New 
scenarios  can  be  created  by  entering  user  chosen  target  parameters  in  the  screen  for  targets  or  by 
accepting  preloaded  scenarios.  The  information  fusion  algorithm  can  be  executed  by  means  of  a 
user-interactive  set  of  questions  and  responses  for  selecting  sensor  and  target  parameters.  The 
codes  consists  of  a  set  of  executable  routines  for  generating  tracks,  performing  track  fusion,  and 
for  graphical  display  of  the  results.  Software  prototype  developed  in  this  contract,  is  portable  to 
any  platform  supported  by  .VIATLAB.  including  PC  and  LFNIX  workstations. 
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Figure  1.  Airborne  Surveillance  Scenario 


6.1  SIMULATION  OF  WEIGHTED  COVARIANCE  FUSION  ALGORITHM 


The  weighted  covariance  KTF  algorithm  described  in  Section  5.2  is  implemented  in 
MATLAB  4.0  under  Windows  3. 1  environment.  This  software  has  a  user-friendly 
interface  (GUI)  for  creating  scenarios  and  .selecting  parameters  for  simulating  sensor 
models.  The  GUI  is  written  in  Visual  Basic.  The  executable  is  MTTF.EXE.  This 
software  expects  the  MATLAB  .system  be  installed  in  the  C;\MATLAB  directory.  It  also 
assumes  that  the  file  MATLABRC.M  is  stored  in  that  directory.  Type: 

cd  c:\matiab\fusedat\demo 

fusion 

Note:  The  directory  c:\matlab\fusedat\demo  is  assumed  to  contain  the  data  files  (input) 
defined  by  the  MTTF  program  as  well  as  the  tracking  and  fusion  m-files.  This  is  also  the 
recommended  location  for  the  MTTF.EXE  program. 


.^0 


In  order  to  run  the  MTTF  program,  follow  the  steps  given  below; 


1  In  the  File  Manager,  open  the  directory  c:\inatlab\fusedat\demo  and  run  the 
MTTF.EXE  program. 

Windows  opens  the  Multi.sensor  Track  to  Track  Fusion  dialog  box  as  shown  in  Figure 

2, 

2  To  accept  all  default  settings,  click  the  OK  button.  All  input  data  will  be  read  from 
pre-defined  files  stored  in  the  working  directory  (e.g.  in  c:\matIab\fusedat\demo). 
Notice  that  MTTF  picks  up  initially  the  directory  in  which  it  is  installed.  If  the 
MATLAB  files  are  located  in  a  different  directory,  select  the  directory  using  the 

Working  Directory  tree  box. 

3  To  create  custom  scenarios  change  sensors  settings  in  the  Sensor  1  and  Sensor  2 
boxes.  In  order  to  redefine  other  sensor  parameters  the  option  Build  has  to  be 
selected.  To  generate  randomly  target  trajectories,  in  the  Trajectories  box,  select  the 
Generate  option. 

Upon  pressing  (clicking)  the  OK  button,  the  program  opens  the  Setup  for 
Trajectories  dialog  box  (if  the  Generate  option  is  selected)  as  shown  in  Figure  3.  or 
it  goes  directly  to  the  Sensor  Type  n  parameter  definition  dialog  box  (if  at  least  one 
Build  option  is  chosen)  as  shown  in  Figure  4. 

The  Setup  for  Trajectories  dialog  box  shown  in  Figure  3.  lets  an  user  to  determine 
the  number  of  targets  (default  =  4),  the  sampling  intervals  (default  =  10),  dimension 
(default  =  3-D),  and  initial  position,  initial  velocity,  and  noise  variance. 

The  Sensor  Type  n  dialog  box  shown  in  Figure  4.  (for  Laser  Radar)  permits  an  user 
to  define  such  sensor  parameters  as  Location,  Accuracy,  Minimum  and  Maximum 
Field  of  View,  SNR,  Threshold,  etc. 
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Upon  completion  of  all  setup  tasks.  MTTF  invokes  the  MATLAB  program.  Through  the 
MATLABRC.M  file,  the  working  directory  is  selected  and  the  FUSION. M  program 
started. 
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Figure  2.  Sensor-Scenario  Definition  Dialog  Box 
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Figure  3.  Trajectory  Setup  Dialog  Box 
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Figure  4.  Laser-Radar  Definition  Dialog  Box 
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Three  routines  were  developed  to  perform  Monte  Curio  simulution  of  two-sensor  tiuck 
fusion.  Names  of  these  routines  are:  demo67,  demo68,  and  demo78.  The.se  routines  can 
be  executed  by  typing,  say  "demo67<CR>".  Details  of  these  routines  are  given  below. 
Monte  Carlo  simulation  was  performed  over  a  time  period  of  100  sec.  Figure  5.  shows 
the  four  targets  and  the  IR  tracks  of  those  targets.  Figure  6.  shows  the  position  errors 
with  this  tracker  from  t  =0.0  to  t  =  100  sec.  Since  the  2-D  IR  sensor  is  tracking  a  3-D 
target,  position  errors  increa.se  with  time.  Figure  7.  shows  the  Radar  tracks  and  the 
position  errors  with  this  3-D  sen.sor  are  shown  in  Figure  8.  In  this  case,  the  position 
errors  are  less  than  the  IR  tracker  errors.  Figure  9.  shows  the  Laser  Radar  tracks  and 
Figure  10.  shows  the  position  errors  for  this  tracker.  As  expected,  position  errors  near 
target  cross-over  are  larger  than  in  rest  of  the  scenario.  Results  of  fusion  of  these  tracks 
using  the  weighted  covariance  KTF  algorithm  are  given  below: 

1 )  demo67  implements  fusion  of  tracks  from  IR  and  Radar  sensors.  Figure  1 1 .  shows  the 
fused  tracks  and  Figure  12.  shows  the  fusion  errors.  Figure  12.  shows  that  considerable 
reduction  in  position  errors  can  be  achieved  by  performing  track  to  track  fusion. 

2)  demo68  implements  fusion  of  tracks  from  IR  and  La.ser  Radar  sensors.  Figure  13. 
shows  the  fused  tracks  and  Figure  14.  shows  the  fusion  errors.  Figure  14.  shows  that 
considerable  reduction  in  position  errors  can  be  achieved  by  performing  track  to  track 
fusion. 

3)  demo78  implements  fusion  of  tracks  from  Radar  and  Laser  Radar  sensors.  Figure  15. 
shows  the  fused  tracks  and  Figure  16.  shows  the  fusion  errors.  Figure  16.  shows  that 
considerable  reduction  in  position  errors  can  be  achieved  by  performing  track  to  track 
fusion. 
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Figure  6.  IR  Tracker  Position  Errors 
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Figure  10.  Laser  Radar  Tracker  Position  Errors 
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and  Laser  Radar 


Figure  16.  Fused  Position  Errors  for  Radar  and  Laser  Radar 


6.2  SIMULATION  OF  INFORMATION  FUSION  ALGORITHM 


Track  fusion  algorithm  employing  information  fusion  methodology  was  implemented  in 
MATLAB  4.0  under  Windows  3. 1  environment.  This  software  can  be  installed  by 
creating  a  directory,  for  example:  c:\matlab\fusedat\nova,  and  by  copying  all  the  files  in 
the  diskette  into  the  directory.  To  start,  double  click  on  the  MATLAB  icon  and  the 
command  window  appears.  Type  "cd  c:\matlab\fusedat\nova"  and  the  program  can  be 
executed  interactively. 

Three  routines  were  developed  to  perform  Monte  Carlo  simulation  of  two-sensor  track 
fusion.  Names  of  these  routines  are:  demo7,  demo8,  and  demo9.  These  routines  can  be 
executed  by  typing,  say  "demo7".  To  use  the  default  parameters  in  these  demos,  hit 
<CR>  on  every  pause/prompt  for  input.  The  default  parameters  are  listed  in  the 
parenthesis  after  each  prompt  for  input.  When  multiple  numbers  are  displayed  in  a  single 
request,  they  are  the  only  allowable  inputs  and  the  last  number  is  the  default  when  no 
input  is  given.  Details  of  the  three  demos  are  given  below: 

1)  demo7  implements  fusion  of  tracks  from  IR  and  Laser  Radar  sensors.  This  demo 
illustrates  the  results  of  fusion  of  a  4-D  sensor  (Laser  Radar)  and  a  2-D  sensor 
(IR).  A  centralized  fusion  architecture  is  used  where  the  tracks  from  the  two 
sensors  are  processed  and  fused  at  a  central  fusion  center. 

2)  demo8  implements  fusion  of  tracks  from  IR  and  Radar  sensors.  This  demo  illustrates 
the  results  of  fusion  of  a  2-D  sensor  (IR)  and  a  3-D  sensor  (Radar).  A  centralized 
fusion  architecture  is  used  where  the  tracks  from  the  two  sensors  are  proces.sed  and 
fused  at  a  central  fusion  center. 

3)  demo9  implements  fusion  of  tracks  from  Radar  and  Laser  Radar  sensors.  This 
demo  illustrates  the  results  of  fusion  of  a  4-D  sensor  (Laser  Radar)  and  a  3-D 
sensor  (Radar).  A  centralized  fusion  architecture  is  used  where  the  tracks  from 
the  two  sensors  are  processed  and  fused  at  a  central  fusion  center. 


41 


63  PERFORMANCE  COMPARISON 


To  compare  the  performance  of  the  two  fusion  algorithms,  the  example  used  in  [3]  is  used  as  a 
basis  for  comparison.  In  this  example,  a  single  target  is  tracked  by  two  sensors.  The  kinematic 
model  of  the  target  is  given  by 


W(tk) 

with  sampling  rate  AT  =  I  sec.  The  input  noise  W(t|J  is  assumed  to  be  zero  mean  white 
Gaussian  whose  variance  is  q. 

Observations  from  the  two  sensors,  denoted  by  the  superscript  m,  are  modeled  as 

Z'"(t,j)  =  [l  0]X(tk)  +  v™(tk)  form=l,2 

where,  sensor  observation  noises  are  assumed  to  be  white  Gaussian  and  independent  of  each 
other.  For  the  sake  of  simplicity,  it  is  also  assumed  that  the  sensors  are  similar  so  that 
Var[v‘(t|J]  =  Var[v^(tk)] »  1.  To  analyze  the  performance  of  the  two  algorithms,  the  fused 
covariance  equation  (5o20)  can  be  written  as 

P^tfc  /t,)-'  =2P^tk  /t,)-'  -2p5(t„  tP^dj  /!,.„) 

where,  for  identical  sensors  =  P'ft^/tij)  =  P'lt^/ti^)  is  the  track  variance,  and 


1  I 
0  1 


X(tk)  + 


1/2 

1 
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are  the  extrapolated  covariances  of  the  individual  trackers  and  fused  tracks  respectively,  and  n  is 
the  number  of  sampling  intervals  between  communications  (fusion).  The  steady  state  fused 

covariance  can  be  obtained  by  letting  P^Ct^  /  ti^)  =  P^(t|c.n  /  t^.n)  =  P*  and 

P*^(tk  /  tfc)  =  P*^(tk.n  /  tk-n)  =  P  and  solving  for  components  of  P. 


Figure  1 7.  shows  the  ratios  of  the  elements  of  the  covariance  matrix  P  and  P*  for  a  wide  range 


of  values  of  the  variance  of  the  process  noise  qr,  where 


pi  p2 
p2  p3 


Figure  18.  shows  the 


ratios  of  the  areas  of  the  ellipses  of  uncertainty,  which  are  proportional  to  the  square  root  of  the 
covariance  matrices.  The  dotted  lines  are  the  results  based  on  (2-10),  and  the  solid  lines  are  the 
results  based  on  (5-20)  for  several  fusion  intervals  denoted  by  n.  These  figures  demonstrate  that 
the  information  fusion  algorithm  (5-20)  is  sensitive  to  the  level  of  process  noise  (particularly  for 
p3)  and  the  length  of  communication  interval.  This  is  understandable,  because  when  the  input 
noise  variance  is  increased  or  when  the  time  interval  between  fusion  is  longer,  the  impact  of  the 
underlying  process  noise  will  be  significant.  Figure  17.  and  Figure  18.  also  show  that  when  the 
process  noise  variance  or  the  communication  interval  is  very  large,  ratios  of  elements  of 
covariance  matrices  approach  0.5.  This  is  equivalent  to  the  case  where  the  cross-correlation 
matrix  P^  =  0. 


It  should  also  be  noted  that,  even  though  the  weighted  covariance  algorithm  (2-9)-(2-10)  is 
inexact,  the  error  is  small  and  is  fairly  insensitive  to  the  level  of  q.  This  approach  to  track 
fusion,  however,  requires  sensors  from  the  remote  stations  to  transmit  at  every  t^  ,  the 
observation  matrix  Flft^)  and  the  Kalman  filter  gain  matrix  Kft^)  in  addition  to  the  track  and  its 
variance.  Since  the  existing  communication  links  have  limited  bandwidth,  it  may  be  useful  to 
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perform  a  trade-off  study  involving  the  improvement  due  to  fusion  and  the  link  bandwidth  for 
the  two  approaches.  Such  a  trade-off  study  could  also  consider  the  possibility  of  transmitting 
raw  measurements,  instead  of  the  processed  tracks,  to  a  fusion  center.  Fusion  of  raw 
measurements  from  the  remote  stations  leads  to  optimal  track-to-track  association  and  fusion 
[ll]. 

When  n=l  (full-rate  communication),  the  information  fusion  algorithm  (5-20)  is  optimal  and  the 
results  are  identical  to  the  one  obtained  by  “measurement  fusion”  [15].  Figure  17.  and 
Figure  18.  provide  analytical  performance  trade-offs  for  designing  a  fusion  system  under  various 
operating  conditions.  For  example,  when  the  process  noise  is  relatively  small,  the  information 
fusion  algorithm  is  a  good  choice  because  it  provides  almost  optimal  performance  with 
minimum  communication  requirements.  On  the  other  hand,  when  the  process  noise  is  large,  the 
weighted  covariance  fusion  approach  (2-9)-(2-10)  is  preferable  because  it  incurs  minor 
performance  degradation. 
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Figure  17.  Ratios  of  Elements  of  Covariance  Matrix  for  Fused  State  Estimates  to 

Single  Sensor  State  Estimates 
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Figure  18.  Ratios  of  Area  of  Error  Ellipses  for  Fused  State  Estimates  to 
Single  Sensor  State  Estimates 
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SECTION  7 


SUMMARY 

This  report  describes  the  development  of  improved  track  to  track  fusion  algorithms  for 
C^I  employing  multiple  dissimilar  sensors.  Three  sensor  models  Radar,  IR,  and  Laser 
Radar)  were  generated  using  Monte  Carlo  simulation.  These  sensors  detect  four  targets, 
perform  report  to  track  correlation  and  create  tracks  at  different  data  rates.  For  this  report 
it  is  assumed  that  the  Radar,  IR,  and  Laser  Radar  sensors  create  tracks  at  rates  4  sec., 

2  sec.,  and  I  sec.,  respectively.  These  tracks  are  synchronized  by  using  the  KTS 
algorithm.  The  synchronized  tracks  are  associated  by  employing  the  KTA  algorithm. 

This  algorithm  associate  track  pairs  by  computing  minimum  closeness  score.  As  a 
measure  of  performance,  probability  distributions  of  correct  track  association  and  false 
track  correlation  are  computed  by  incorporating  the  cross-covariance  in  the  test  statistic 
for  association.  A  trade-off  study  involving  the  size  of  the  association  gate,  probability  of 
correct  association  and  probability  of  false  correlation  was  also  performed.  Associated 
tracks  are  then  fused  using  the  KTF  algorithm.  This  algorithm  weights  the  tracks  by 
certain  covariances  which  are  functions  of  track  covariances  as  well  as  their  cross¬ 
covariance,  and  then  summed.  Numerical  simulation  shows  that  the  position  errors  of  the 
individual  trackers  can  be  significantly  reduced  by  performing  kinematic  track  fusion 
which  results  in  more  efficient  and  improved  operation  of  airborne  surveillance  systems. 

In  addition,  a  RTF  algorithm  based  on  information  fusion  approach  was  also  developed. 
This  algorithm  is  efficient  because  its  implementation  requires  less  computation  but  does 
not  incorporate  the  cross-covariance  between  the  candidate  tracks  for  fusion. 

Performances  of  the  two  algorithms  are  compared  by  means  of  covariance  analysis. 
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APPENDIX  A 


DESCRIPTION  OF  SENSOR  MODELS 

This  section  describes  the  panuneters  of  the  three  sensors  (Radar,  IR,  and  Laser  Radar) 
which  are  used  to  detect  the  targets.  In  this  report,  these  sensors  are  assumed  to  be 
colocated.  After  detection,  the  sensors  perform  report  to  track  correlation.  These  tracks 
are  transmitted  to  a  central  fusion  center  where  track  fusion  is  performed. 

Radar  Sensor  Model 

The  radar  model  used  in  this  simulation  is  described  in  [12].  Output  of  the  radar  is  the 
observation  vector  [Range,  Azimuth,  Elevation].  The  equation  describing  the  range  and 
signal  to  noise  ratio  is  given  by 


R 


PAqL 

16(kT)(S/N) 


(t./Q)'"* 


Nominal  values  of  the  critical  parameters  are  used  in  this  study  are  given  below: 


P  -  500  Watts 

D  (diameter  of  aperture) »  Ift 
A  » 7t(l)^/  4(3.28)^]  »  0.073  m^  aperture  area 
k  =  1 .38x10'^  Joule/’K  *  Boltzman  constant,  T  *  600®K 
<T  =  1 0  m*  ( small  aircraft) 

L  =  L^3,  X  L,u„os.  “  Radar  (transmit/receive)Losses  =  0. 1 
Latinos  “  Atmospheric  Losses  (from  3GHz  to  8GHz)  =  0.6 
L  *0.1  X  0.6  *0.06 

Q  =  Angular  scan  coverage  *  100  x  (120®)  x  (it/l80)^  *  0.365  steradians 
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tj  =  Scan  duration  =  3  sec. 

X(at  8  GHz)  =3xl0*/8xl0’  =  0.0375  m 

0BW  =  0.0375/0.305  =  0. 123  radians  =  1^  (elevation) 

Assuming  a  circular  aperture,  null-null  beamwidth  =  2x7®  »14® 
9bw  (azimuth)  =  2° 


IR  Sensor  model: 

The  IR  sensor  model  used  in  this  simulation  is  described  in  [13].  Output  of  the  IR  sensor 
is  the  observation  vector  (Azimuth,  Elevation).  The  equation  describing  the  range  and 
signal  to  noise  ratio  is  given  by 

r  TtPJP  L  ,i/2.at3vi/4 

Uim)iS/W 

Nominal  values  of  the  critical  parameters  are  used  in  this  study  are  given  below 

D  =  Optical  aperture  =  1  ft.  =  30.48  cm. 

D*  =  Dielectric  constant  =  2X10’®cm-HZ''^“/Watts 
L  =  System  and  atmospheric  losses  *  0.1, 

n  =  Number  of  resolution  elementsfdetectors)®  3600  (600  in  elevation  and  6  in  azimuth) 
f#  *  Focal  number  *  1.5,  t,  =  Frame  time  ®  2  sec. 

Q  =  Total  solid  angle  coverage  =  120®(ele)xl0®(azi)x(7t/180)^  =  0.365  steradians 

J  =  Source  radiation  intensity  =  40  watts/  steradians 

©Bw  ~  Beamwidth  =  7JD  =  10'*  m/  0.3048  m  =  32.8  m  rad  =  0.002  degrees 

Resolution  in  elevation  =  angular  coverage/  number  of  detectors  =  10®/600  =  0.0166  deg 

=  0.29  m  rad 
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Laser  Radar  Sensor  Model 


The  Laser  Radar  sensor  model  used  in  this  simulation  is  described  in  [14].  Output  of  the 
Laser  Radar  sensor  is  the  observation  vector  (Range,  Range  Rate,  Azimuth,  Elevation). 
Sensor  model  given  below  is  based  on  the  assumption  of  a  CO2  laser.  For  a  coherent 
receiver,  the  equation  describing  the  range  and  signal  to  noise  ratio  (for  a  point  target)  is 
given  by 


r*  <TnSYS’lATM  J/4 

Nominal  values  of  the  critical  parameters  are  used  in  this  study  are  given  below 

R  ®  system  range  to  target  *  10’  meter 
Et  =  transmitted  laser  energy  »  I 
S/N  *  signal  to  noise  ratio  =  3 1 .4 
D  =  receiver  aperture  diameter  =  O.lm, 
riATM  =  atmospheric  transmission  factor  =  0.1 
TisYS  “  system  transmission  factor  =  0. 1 
N  ^  Quantum  noise  of  coherent  optical  receiver 
h  •  Planck's  constant »  6.626x10*^  J-s 
X  ■  Wavelength  of  radiation  *l0nm 

.  _  .  CT 

Standard  deviation  of  the  measurement  error  in  Range  is  =  — ————— 

2^2(5 /N) 


Standard  deviation  of  the  measurement  error  in  Doppler  is 


V3X 

*7rcV2(S/N) 
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Standard  deviation  of  Azimuth  and  Elevation  is  = 

where,  t  =  pulse  width.  For  example,  the  nominal  system  parameters  of  S/N  =  50, 

T  =  10"^  sec.,  A.  =  10'^  meter,  and  D  =  O.l  meter  would  give  rise  to  the  following 
measurement  noise  standard  deviations:<TR  «  15  m,  ao  «  0.5  m,  ©a  ~  ^ 

For  the  targets  within  the  sensor  field  of  view,  the  detection  model  [14]  can  be  described 
as 

_j _ 

p  _pl+{S/N) 

where.  PpA  is  the  probability  of  false  alarm. 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Sun/eillance, 
Communications,  Command  and  Control,  Intelligence.  Signal  Processing, 
Computer  Science  and  Technology,  Bectromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


